Low-bias mobility vs high-bias saturation current in graphene based transistors 
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We describe the fundamental trade-offs in engineering the mobility, current saturation and ON- 
OFF ratios in graphene transistors. Surprisingly, the trade-offs arise solely from an asymptotic 
constraint on the high energy bandstructure and independent of scattering processes. This places 
graphite derivatives (bulk monolayer graphene, uniaxially strained graphene nanoribbons, carbon 
nanotubes and bilayer graphene) on the same 3-parameter mobility-bandgap-scattering length — 
Egap — X) plot, proximal to other semiconductors. In addition to this low-bias trade-off, the high bias 
current bears signatures of the underlying saturation mechanism, arising through phonon scattering 
or F-point suppressed density of states opening bandgap. 



The incredible properties of graphene have stimulated 
intense exploration into its potential as an electronic de- 
vice [ij. Graphene retains many admirable properties of 
nanotubes, such as a high mean-free path for scattering 
and a paucity of surface dangling bonds, minus many of 
their disadvantages, such as insensitivity to atomic chi- 
rality and the possibility of top-down lithographic pat- 
terning. This has led to the study of graphene devices 
that range from evolutionary (such as monolithically pat- 
terned wide- narrow- wide ribbons to revolutionary 
(such as Veselago 'waveguides' and switches^ , and Bi- 
layer pseudospin field effect transistors or BisFETs 

The biggest problem for a graphene-based electronic 
switch is undoubtedly its narrow band-gap. Various ef- 
forts to engineer the band-gap have consistently seen gaps 
narrower than ^ 200 meV. Such a small band-gap threat- 
ens to severely limit the room temperature ON-OFF ra- 
tio to only ~ 50 (keeping in mind that for switching 
applications the drain and gate bias voltages are typi- 
cally comparable). The small bandgap also generates a 
large subthreshold swing ~ 145 mV/decade |5| and large 
source-to-drain tunneling currents at modest voltages [6| . 
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The reliance on chemical means or quantization for 
band-gap creation compromises overall device scalability. 
More seriously, the bandstructure of graphitic derivatives 
pose a stringent restriction on their electronic properties, 
so that opening a band-gap via strain or chemistry (as 
in graphane [7|) undermines the incredible mobility of 
graphite electrons arising from their small effective mass. 

In this paper we explore fundamental band-related con- 
straints on both the low-bias mobility and the high- 
bias current saturation in FETs created out of graphite- 
based materials. We outline an intimate trade-off be- 
tween mobility (switching speed) and ON-OFF ratio 
(bit error rate) that stems from an asymptotic con- 
straint on the high-energy bandstructure of all the 
well known graphite derivatives, including single layer 
graphene (SLG), strained graphene nanoribbon (sGNR), 
carbon nanotubes (CNT) and bilayer graphene (BLG), 
even without a corresponding change in scattering length. 
Lastly, we outline different mechanisms behind current 
saturation in graphene, prior to band-to-band tunneling - 



FIG. 1. Regardless of the bandgap opening mechanism, all 
monolayer graphitic materials such as SLG, CNTs, GNRs 
are constrained to a single (a) effective mass versus bandgap 
curve. Graphitic materials are fundamentally limited by a (b) 
dispersion relation that is pinned to its linear form at high 
energies. Compared to other common group IV and III-V 
semiconductors at corresponding energy bandgaps, graphene 
electron carriers are not much lighter. 

focusing on discernible inelastic electron tunneling spec- 
troscopy(IETS) Q signatures arising from saturation due 
to the F-point, bandgaps and phonon scattering. 

Relating mobility to bandgap. Irrespective of the un- 
derlying bandgap opening mechanism, the bandstruc- 
ture of all graphite derivatives reach a linear dis- 
persion at higher energies as expressed in E{k) « 
±y/{Egap/2y + {hvoky, where h is reduced Planck's 
constant. For SLG, vq = 3aot/2h is the velocity of the 
high energy graphene electrons 10* cm/s, with ap be- 
ing the C-C bond-length and t ^ 2.5 eV being the bond 
energy. Whether this slope itself changes depends on the 
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wavelength of the perturbing potential (e.g. confinement 
or strain) relative to the graphite lattice constant. More 
specifically, 
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Meanwhile the terms for sGNR [9| dispersion, A', B' C', 
7i, 73, p, and (N+1) are defined in the reference [10| . 
Setting appropriate strain terms to zero recovers the dis- 
persion for GNRs with first and third nearest neighbor 
interactions, and edge distortion For BLG the terms 
A, B, C, and a are defined in reference [l2j. Since the 
opening of Egap does not influence its high energy sec- 
tor that is pinned to the linear dispersion, it flattens the 
band- curvature and increases its efl^ective mass (Fig. 1). 
There is thus an inherent trade-off between ON-OFF 
ratio (bandgap) and mobility (effective mass) based on 
bandstructure considerations alone. Taking the second 
derivative at k = gives the effective mass at the band- 



bottom, TOq 



Egap/vQ, implying that the kinetic energy 



gained by the electrons and holes equals the energy lost 
from the crystal potential during band-gap opening. The 
definition of effective mass away from the band-bottom 
requires careful consideration. One can write the mobil- 
ity /X = q\/m*v, where the effective mass to* and the 
carrier velocity v are in general, energy- dependent, and 
A is a constant mean free path. Invoking the dynamical 
definition of effective mass, /x can be rewritten in terms 
of the gate voltage by first replacing to* v with hk p at the 
Fermi energy (i.e., using to* = p/v reduces to TOq only 
at the band-bottom [l3|). From the graphene dispersion, 
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for the electronic sector. 
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where Ep 
We then have 
qXvn 



Egap/2+qVG 



V^VcEg 



(qVa) 
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Since Vg also determines the 2-D electron density 
through kp = y^m2zJ, the gate dependence of the mo- 
bility translates to a dependence on n2D ~ Vg{Vg + 
Egap) I L-kH^Vq for a channel of length L. 



While much interest in graphene, such as for RF ap- 
plications, stems from its ultra-high mobility [l^ . for dig- 
ital switching applications one gets a more complete pic- 
ture by plotting an Egap vs curve for various As, as 
in Fig. 2. One can see that for zero band-gap at room 
temperature, we get the maximum theoretical mobil- 
ity qXvo/kT ~ 400,000 cm^/Vs for a 1.2/xm scattering 
length. The perceived advantages of graphene get com- 
promised when one factors in both its switching speed 
and its ON-OFF ratio. Scattering processes from impu- 
rities [H, [3l , bulk or interfacial phono ns [l^ E^l : and 
line edge roughness [l^ 21, 22 1 in GNRs are expected 



to shorten A, further decreasing the mobility in addition 
to the fundamental bandstructure constraints. From our 
analysis it seems that once the advantages of low m* and 
high A are compromised, graphitic switches act compa- 
rable if not worse than Indium Antimonide (InSb) and 
Gallium Arsenide (GaAs), and marginally better than 
Silicon (Si). Using the dispersions (1) and the generic 
(/i — Egap — A) formulation (2), we show universality 
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FIG. 2. 3- parameter (jj, — Egap — A) illustrates the fun- 
damental cost in opening a bandgap to carrier mobility in 
graphitic materials. With an electron density of 2xl0^^c?n^'^ 
and A « 1.2fj,m, zero-bandgap graphene reaches a mobility of 
230, OOOcm^/ys matching experiments [3. Compromising A 
as mentioned in the text potentially lowers carrier mobility 
comparable to some group IV and III-V semiconductors [la] 
for a given Egap- (b) All monolayer graphitic derivatives fall in 
a single curve with small difference in the value of vq [ii|]. For 
the same set of armchair (3p+l, 0) GNRs, we show in the sub- 
plot that uniaxial strain((T — 5%) can close bandgaps(extend 
fi). BLG trend deviates from monolayer graphitic derivatives 
due to its non-linear zero-bandgap dispersion and its unique 
band-bottom features in the presence of a bandgap [l^ . [l7| . 
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in trade-offs amongst single layer graphitic derivatives 
(Fig. 2b). Wc would emphasize the trade-offs in open- 
ing Egap on fi is synonymous with conductance((j' = 
^M{E)Tj^), where the mode density (M(E)) near 
hand-bottom is diluted independent of A, while transmis- 
sion per mode(T) and graphene length (L) are constants. 

High bias current: Modeling scattering. We saw that 
the low-bias mobilities of all graphitic derivatives are 
intimately connected with their band-gaps. A similar 
universality arises for the high-bias current. The dimin- 
ishing density of states for graphene near the F point 
conspires together with various scattering mechanisms to 
create a current-voltage (FV) characteristic that shows 
a tendency to saturate, followed sharply by a rise due to 
band-to-band tunneling. A closer look at the FV and its 
derivatives provides useful insights into the underlying 
scattering mechanism. To this end, we will now discuss 
how saturation occurs in our model for graphene FETs. 

Electron transport through a nanoscale object is 
described using the non-equilibrium Green's function 
(NEGF) formalism that simplifies to Fandauer theory 
for coherent quantum flow [23|, 
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where /i_2 are the bias-separated contact Fermi-Dirac 
distributions, the transmission T ~ irace(FiGF2G^), 
with the Green's function G = {ES - H - S, H 

and S being the channel overlap matrix, Hamiltonian and 
self-energy matrices respectively, and F.^ = i(S]j — s|) is 
the level broadening. A recursive technique allows us 
to get E for any layered contact structure in the non- 
interacting limit. The influence of scattering modifies 
the Green's function through Dyson's equation: 
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where Go is the unperturbed response and Eg repre- 
sents perturbing interactions from scattering processes 
such as from impurities (elastic) or phonons (inelas- 
tic). From an energy dependent density of states 
Ec) + Q{E^ — E)] we can extract 



DiE) = ^mE 



the imaginary part of Go and Hilbert transform to get the 
real part as well. S is the area of the SLG sheet (500nm 
X 500nm). For inelastic phonon scattering, we can invoke 
the self-consistent Born approximation as follows: 
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where E™2°"*(-^) ^'^^ the self-energy terms for the con- 
tacts and Ep'^'°^* is the self-energy term for the inscat- 
tering and outscattering that occur at a specific en- 
ergy (E=:±?iaj). The scattering deformation potential, 
'DQ{uj)—Xi,pcfiuj [23 |. where \epc is the electron-phonon 
coupling constant (epc). For acoustic phonons we use 
a Debye model to cover the range of frequencies up to 



^huj^/kBT 



is the 



the Dcbey frequency. Na{uj) 
phonon Bose-Einstein distribution describing the popu- 
lation of phonons at hu, and /j is the terminal current. 

Current Saturation Mechanisms. Prior to band-to- 
band tunneling in a gFET is an onset of scattering- 
induced-saturation (SIS) due to phonons from lattice vi- 
brations either inherent to graphene or the underlying 
substrate atoms. We focus on intra-valley longitudinal 
optical (FO) and longitudinal acoustic (FA) phonons as 
examples that strongly couple with electrons and infiu- 
ence transport. Using the Einstein model FO phonons, 
are active at single energy huj w 0.19eV, while FA 
phonon are approximated with the Debye model are ac- 
tive across a bandwidth of energies {hoj « OeV^ — 0.16eF). 

Optical phonons are known to saturate the current in 
metallic CNTs depending upon on the strength of Vojuj) 
Hi]. Va{Lu) for both FA and FO phonons from CNTs Q 
are applied to our scattering model for SFG without loss 
of generality to the onset of SIS on current. Accordingly, 
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FIG. 3. IV-characteristics and lETS for scattering induced 
saturation (SIS) in zero-bandgap gFETs (red-curves). Bal- 
listic case in dashed-blue curves, (a) LO phonons in the IV 
(Vg = 1.5V, 1.8V) are active at &j f» 0.19eV. (b) Vibronic 
signatures in lETS at 0.19eV are seen with overtones at higher 
harmonics, (c) LA phonons are active over a range of ener- 
gies from huu « OeV — O.lGeV and contribute to decrease in 
conductance (Vg = 1.25V, 1.50V, 1.75V) (d) Corresponding 
lETS shows a broadened valley near the Debye frequency for 
LA modes. *Parameters : F, ~ 2eV [2^, graphene area 
(S): 500nm x 500nm, LO phonon Aepc ~ 0.804, LA phonon 
Aepc ~ 0.40 . Note different voltage ranges in lETS (b,d). 
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FIG. 4. Narrow graphene bandgaps difficult to discern in the 
IV create a well defined plateau in the lETS. Zero-bandgap 
gFET in absence of SIS has a linear lETS fiom the symmet- 
ric infiection in the IV at the F-point. The diamond sym- 
bols mark the F-point shift at different levels of electrostatic 
doping(Ss„p = 0.4eV^, Vg = 1.5\/; E.^r. = 1.0\/, Vq = l.QV ). 

the onset deviation from the balhstic-IV is perceived at 
energies where phonons are activated (Fig. 3 ). 

We thus find that a critical feature in LO phonon in- 
duced saturation is that the onset voltage for saturation is 
pinned to a gate-independent value given by the phonon 
frequency. This feature can be accentuated by measur- 
ing the second harmonic response and plotting the cur- 
vature of the I-V curve, as in lETS measurements. A 
sharp valley in the d^I/dV^ curve is expected to arise 
at a gate independent voltage for LO phonon SIS (Fig. 
3b). In addition, higher order harmonic side-bands in 
the lETS (Fig. 3b) are a visible and reminiscent feature 
in CNTs [2J|. Meanwhile, the strength of LA phonons 
increases with D{E) and are averaged over a range of fre- 
quencies creating a broadened valley in the lETS whose 
minimum is short of the Debye frequency (Fig. 3d). 

However, we add that opposite effects of scattering 
is visible in weakly coupled contact-channel gFET. Low 
bias regimes with a paucity of states near the T-point are 
prone to forward scattering which enhances conduction 
and further evidenced by peaks in the lETS at the rele- 
vant vibronic frequency. The opposite effect prevails at 
higher bias regimes with more states. 

For zero-bandgap graphene, the paucity of states near 
the F-point creates a seemingly short plateau that is re- 
ally an inflection in the IV. We denote the linearly de- 
creasing D{E) near the F-point at low bias regimes as the 
source of F-point induced saturation (GIS). To achieve 
currents independent of voltage requires a bandgap in 
graphene. However bandgap induced saturation (BIS) 
for even narrow bandgap graphene is nearly indiscernible 
from the IV alone due to band-to-band tunneling. Al- 
though difficult to measure, lETS should still resolve a 
short plateau corresponding to a narrow bandgap (Fig. 4) . 

Conclusion. In this work we point out that attempts 
to engineer a bandgap in gFETs invariably reduce their 
mobilities, based on asymptotic band-related constraints 
alone. On one hand, gFET as a switch is fundamentally 
compromised by ^—Egap—X relation. On the other hand, 
subthreshold conduction from band-to-band tunneling in 
narrow -bandgap (higher-mobility) gFETs compromises 



ON-OFF current ratios, thus requiring a careful balance 
in future design considerations. Finally, we differentiate 
IV saturation in gFETs arising from bandgaps, F-point 
and vibronic scattering. Based on these studies, we be- 
lieve it should be straightforward to create a compact 
model for IVs based on any gFETs, that can incorporate 
a multitude of geometrical and atomistic variations as 
well as underlying physical scattering mechanisms. 
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